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Nuclear Thermal Rocket Introduction 


Few rocket propulsion concepts offer the combination of high thrust and reasonable efficiency that can be obtained 
from a Nuclear Thermal Rocket (NTR). Long considered one of the most basic forms of Nuclear Propulsion, the solid-core 
nuclear thermal rocket engine concept typically employs a uranium fueled nuclear reactor core and hydrogen (H 2 ) gas working 
fluid. The H 2 gas acts first as fuel rod coolant as it passes through the nuclear reactor core followed by rocket working fluid 
when the then super heated hydrogen is expanded out of a nozzle in order to produce thrust. 



Famously the Nuclear Engine for Rocket Vehicle Application, or NERVA test program of the 1960's extensively 
demonstrated the feasibility of the nuclear thermal rocket design concept, testing nearly 30 nuclear thermal rocket engines such 
at the KIWI-B4 which produced approximately 75,000 pounds of thrust with an efficiency (isp) of 825 seconds. 


Atypical Nuclear Thermal Propulsion (NTP) 
rocket engine is the product of a host of design 
tradeoffs including but not limited to heat transfer, 
structural including fuel thermal swelling, core 
reactivity, radiation shielding, thrust, volume, 
mass, and etc 

Fuel surface area plays a key role in heat 
transfer along with the availability of coolant 
pathways thus the surface to volume ratio of a 
NTP concept is often used as a quality measure 
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NERVA/Rover NTP limitations 


* NERVA/ROVER designs were very large and very heavy 

* A typical NERVA program Nuclear Thermal Propulsion (NTP) system averaged a thrust to weight ratio of about 4:1 

* Fails to meet the low mass requirements, and limited available space of today's space launch vehicles 

* Precludes most unmanned missions 

* Modern chemical liquid oxygen - liquid hydrogen (LOX/LH2) engines routinely achieve 50:1 and higher 

* Space Shuttle Main Engine T/W is about 65:1 in vacuum 





Review of Modern NTP Concepts 


Modern NTP designs would like to achieve an efficiency (isp) of 1000 seconds versus NERVA 825 seconds 

* NTP fuel rods will be exposed to extreme operating temperatures and pressures to achieve an isp of 1000 seconds 

* Operating in excess of 5000 °f, and 1000 psi 

Next generation designs would like to achieve a superior thrust to weight ratio over NERVA 4 to 1 

* Will require non-traditional fuel rod materials, geometry, and manufacturing techniques 
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Fuel Particle Evolution 


A literature review of several modern solid core potential NTP concepts was performed 

The purpose of the review was to identify both the strong and weak points of each fuel element concept, hoping to glean 
the ideal fuel rod characteristics for optimum thermal heat transfer to the coolant/working fluid while still maintaining 
structural integrity 

* A large fuel rod surface area ideally shaped to promote heat transfer to the hydrogen working fluid while 
simultaneously avoiding structural degradation due to the high operating temperatures and pressures 

* Pebble Bed Reactor 

* Particle Bed Reactor 


• Foam Core Reactor 
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Localized Thermal Instability 

Deleterious thermal effects systemic to Pebble Bed and Particle Bed designs were identified 

* Both designs suffered from unexplained localized thermal instabilities leading to thermal runaway and fuel melting 


Similarities of the two designs 

* Small localized hot spots can occur within a reactor core for any one of a host of reasons unrelated to coolant flow 

* Design flaw 

* Manufacturing flaw 

* Fuel migration within the pellet/particles affecting fission reactivity 

* Small instabilities that are randomly inherent to a large scale fission reaction 

* Coolant flow will resist localized high temperature regions 

* Coolant fluid properties change with temperature 

* Both Pebble Bed and Particle Bed reactors by design utilize randomized coolant channel pathways 

* Traditional designs such as the close packed hexagonal reactor cores have structured coolant flow pathways 


* Have a greater thermal stability 
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Pebble bed Reactor 


Helium cooled, graphite moderated, high temperature reactor 

Use thousands of marble to softball sized ceramic coated uranium fuel pebbles 


Graphite shell Coaled particle 

(Fissile material a. fertile material) 


U0a + ThOj 
or 

UC+ThC 
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The encapsulating style of the pebbles can lead to both 
structural and reactivity issues when built up fission 
product gasses trapped within the pebble by the outer 
coating layers force the uranium fuel to migrate away 
from the center of the pebble 
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Pebble Beds have high operating temperatures (1700 °F helium exit temperatures) compared to most nuclear electric plants 
Pebbles ceramic encapsulating style significantly reduces the negative effects of fuel melting 

The pebbles are stacked in the reactor core like cannon balls resulting in significant volume between them, massively 
increasing the total volume requirements of this type of design and limiting it's applicability to Nuclear Thermal Propulsion 

A 5.1xl0 7 BTU/hr (15 MW) electric, helium cooled pebble bed plant was successfully demonstrated in Germany 

* Suffered from temperature instabilities which resulted in localized thermal runaway for indeterminate reasons 

* Subsequent fuel melting heavily contaminated the pressure vessel with Cs-137 and Sr-90 
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Particle Bed Reactor 


Hydrogen cooled high temperature reactor 


uses billions of 400 micron sized ceramic coated uranium fuel particles held inside a porous walled (Frit) fuel rod 


Fuel Particle Evolution 
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Fuel particles were known to block the porous 
passages of the frit, inhibiting the flow of hydrogen 
coolant 

Reactor operation would result in permanent 
structural deformation (i.e. bowing of the specially 
manufactured porous material frits due to thermal 
expansion) 



* Has an expected efficiency (isp) of 1000 seconds with accompanying operating temperatures of 5000 °F 

* Expected thrust to weight ratio of 35 to 1 

* Testing reveled that marble sized regions spread throughout the fuel rod would suffer from excessive heating, 
subsequent thermal runaway, and localized fuel melting for indeterminate reasons 

* Potentially related to clogging or deformation of the Frit, which would effect reactor criticality 


Foam Core Reactor 


-carbide or carbonitride fuels chemical vapor deposited onto 


Fuel Element comprised of Uranium based bi-carbide, tri 
porous carbon foam 



Excerpt from the original patent application (15) 

The Foam Fuel element eliminates the need for a 
porous (frit) container which are known to 
permanently deformation 

The Foam Fuel element does not have any loose fuel 
particles which are known to block the porous 
coolant passages of other designs 



FIG. illustrates a schematic isometric view of an 
example of a nuclear fuel element 40 comprising porous 
nuclear fuel 44 encased within metal cladding 42, according 
to the present invention. Gas coolant 46 flows through the 
gas-permeable, porous nuclear fuel 44, exchanging heat with 
a high heat transfer efficiency from the high-porosity nuclear 
fuel, due to the large extended surface area of the porous fuel, 
and at a high temperature due to the thimiess of the nuclear 
fiiel itself. 

Excerpt from the original patent application (15) 


Has an expected efficiency (isp) of 925 seconds with accompanying high operating temperatures 
Expected thrust to weight ratio of 35 to 1 



Localized Thermal Instability hypothesis 

A probable cause for the localized thermal instabilities, thermal runaway, and fuel melting was derived 

* Pebble Bed and Particle Bed reactor designs utilize randomized coolant channel pathways that when combined with 
localized high temperature regions would work together to negatively exploit the coolants temperature dependant fluid 
properties diverting it away from where it is needed the most, to cooler, less resistive pathways where it is needed the 
least 

* Given the choice via many optional flow paths the reactor coolant will take the path of least resistance 

* Hot zones offer the highest resistance 

* These minor thermal instabilities will grow over time leading to thermal runaway 

* Thermal Runaways circular effect leads to fuel melting 

* Coolant mass flow rate drops as pellet temperatures increase 

* Pellet temperatures increase as coolant mass flow rates drop 
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Thermal Effects Demonstrator Analysis 

A thermal fluid analysis was performed in support of the NTP concept review in order to verify the root cause of the 
thermal instabilities found in Pebble Bed and Particle Bed reactor cores 

* FLUENT: Thermal Fluid Analysis program 

* PILGRIM : Air Force computational model (grid) generator 


* Thermal fluid analysis of the Thermal Effects Demonstrator 

* Modeling and simulation of the two dimensional "Thermal Effects Demonstrator" model 

* Rapid designing via Pilgrim and FLUENT during trade studies 

* PILGRIM scripting capability allowed for easy grid modification via text editing of the script 

* FLUENT'S capability to couple the analysis of thermal fluid and structural thermal conductivity 

* (Analysis of heat transfer from the fluid "through the wall") 

* Determine the mass flow rate-operating temperature curve of various fuel particles hydrogen exhaust gas 

Limiting Design Factor 

Canned solution designed to test a certain hypothesis as 
simply as possible 

• Time constraints dictated the complexity of the 
model 

• I am not a computational fluid dynamicist, I am 
also not a reactor core designer. I am a thermal 
analyst that identified a potential thermo 
dynamic solution for an as yet unanswered series 
of thermal instabilities that have limited the 
operation of reactors in the past. Test results of 
the hypothesis have revealed surprising answers 
that both computational fluid dynamicists and 
reactor core designers should be made aware of 

• Fuel thermal swelling was not modeled 

• Criticality was not modeled 

• Thermal-Structural-Criticality interactions were 
not modeled (purely a thermal fluid model) 



3" -I- 1 
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Simulation Model 


Physics Model of the Thermal Effects Demonstrator 

* Single velocity inlet boundary condition for H 2 gas 

* Multiple pressure outlet boundary conditions for H 2 gas 

* Volumetric heat generation rate in pellet material 

* Pellets comprised of Molybdenum-Rhenium 

* Heat conduction across particle materials 

* 2-dimensional, Planar analysis 

* Piecewise linear temperature dependant properties (Specific Heat, Thermal conductivity, and viscosity) 


General problem setup information 


* Solver Type 

* Velocity Formulation 

* Time 

* Space 

* Turbulence Model 


Pressure Based 
Absolute 

Steady State Solution 
2-dimensional, Planar 

Spalart-Allmaras (1 eqn) or K-epsilon (2 eqn) with wall functions 


Solution Method 


• Scheme 


SIMPLE Scheme with Pressure-Velocity Coupling 


Spatial Discretization 
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• ModTurbVis 
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Convergence Criteria 

* Low Residuals (<lxl0~ 4 ) with little fluctuation in residual value between iterations 



Computational Model 


PILGRIM 

United States Air Force computational model construction program, or grid builder 
Plot-3d Formatted List (P3dF) 

2-Dimensional Grids 

28,336 Cells ^2 Velocity Inlet 

0.01" wall spacing 
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Case 22 


Case 22 Thermal Fluid Results 


Fixed Temperature assigned to each pellet 

• Pellets assigned a fixed temperature of 2000 °F in run number 1 

9 Center pellet core temperature was incrementally increased in a stepwise manner 

• Increased center pellet temperature in 500 degree increments over 7 runs 

• Reached a center pellet temperature of 5000 °F in run number 7 




Case 22. Run 5 (Center Pellet 4000 °F) 




Case 22, Run 7 (Center Pellet 5000 °F) 
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Case 22 Thermal Fluid Results Cont. 


Case 22 

* Fixed Temperature assigned to each pellet 

* Contour plots show marked increase of coolant velocity in outlier channels as center pellet temperature increases 
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Case 22 Thermal Fluid Results Cont. 


Case 22 

* Fixed Temperature assigned to each pellet 

* Coolant mass flow rate past the overheating pellet is reduced with increasing temperature 


Mass Flow Rate vs Temperature 
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Case 23 


Case 23 Thermal Fluid Results 


Volumetric heat generation rate assigned to each pellet 

• Pellets reached a temperature of 2000 °F (±10 °F) in run number 1 

9 Heat generation rate was increased in the center pellet to increase pellet core temperature 
• Adjusted center pellet temperature in 500 (±10 °F) degree increments over 7 runs 

• Reached a center pellet temperature of 5000 °F (±10 °F) in run number 7 



Case 23. Run 5 (Center Pellet 4000 °F) Case 23. Run 7 (Center Pellet 5000 °F) 
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Case 23 Thermal Fluid Results Cont. 


Case 23 

• Volumetric heat generation rate assigned to each pellet 

• Contour plots show marked increase of coolant velocity in outlier channels as center pellet temperature increases 

• Increased coupling of coolant velocity profiles between neighboring channels is demonstrated 
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Case 23. Run 1 (Center Pellet 2000 °F) 

Case 23. Run 3 (Center Pellet 3000 °F) 
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Case 23. Run 5 (Center Pellet 4000 °F) Case 23. Run 7 (Center Pellet 5000 °F) 
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Case 23 Thermal Fluid Results Cont. 


Case 23 

* Volumetric heat generation rate assigned to each pellet 

* Coolant mass flow rate past the overheating pellet is further reduced with increasing temperature due to thermal 
coupling 

* Analysis demonstrated the boot strap effect, where neighboring pellets which are operating correctly will 
thermally couple with an overheating pellet to create an expanding region of localized thermal instability 


Mass Flow Rate vs Temperature 




Temperature (F) 


18 


Mass Flow Rate (Ibm/s) 


Temperature Dependant Fluid Properties 

Assigning a constant density to the H 2 coolant significantly reduces the negative effect of the thermal instabilities 
Assigning a constant viscosity to the H 2 coolant also reduces the negative effect of the thermal instabilities 

These fixed fluid property results demonstrate that the coolant flow resists localized high temperature regions because of 
changes in fluid properties caused by increases in temperature 

Effect of Const. Density & Viscosity H 2 Coolant on Mass Flow Rate vs. Temp. 
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Results/Observations 


The Thermal Effects Demonstrator does not just mimic the symptoms it also shows why randomized flow paths are vulnerable 
to minor thermal instabilities 

* The difference in pressure along the length of the Thermal Effects Demonstrator (AP > 0) ensures coolant mass flow 

* There is no pressure difference perpendicular to the length (AP = 0) of the main channel 

* When the flow encounters surface area that is perpendicular to the AP at the tips of the pellets it stagnates 

* With no AP to drive the stagnated flow the coolant is instead directed by what begins as relatively small thermal effects 

* These thermal effects then go on to rapidly build up out of control as a pellet is literally starved for coolant 


Observations 



* In an effort to overcome the shortcomings of current NTP designs a great deal of emphasis has been placed on the 
Surface Area to Volume ratio (S/V) of the next generation concepts 

* This has lead to several advanced design concepts that would or could utilize randomized flow paths 

* This analysis has demonstrated that nuclear thermal rocket reactor core designs which utilize randomized flow paths are 
highly vulnerable to the types of minor thermal instabilities that lead to thermal runaway 

* Designs with structured flow paths easily overcome these very same types of thermal instabilities with no impediment 

* A large Surface Area that is Parallel to the AP to Volume ratio (=S/V) positively impacts the design 

* A large Surface Area that is Perpendicular to the AP to Volume ratio (+S/V) negatively impacts the design 

* Parallel Surface Area plus Perpendicular Surface Area equals Total Surface Area (=S) + (+S) = S T 

* A large Total Surface Area to Volume ratio (S T /V) is indeterminate in predicting the thermal stability of a given design 

* A set of cross flow pumps which establish a small AP in the perpendicular direction would easily overcome the initial 
minor thermal instabilities before they have time to build up 

* A cross flow AP may be difficult to apply to highly randomized flow paths of irregularly shaped pellets such as the Particle 
Bed, but should be applicable to Pebble Bed designs 
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Conclusions/Recommendations 

Reactor fuel rod surface area that is perpendicular to coolant flow direction (+S) i.e. perpendicular to the AP creates areas 
of coolant stagnation leading to increased coolant temperatures resulting in localized changes in fluid properties 

* Changes in coolant fluid properties caused by minor increases in temperature lead to localized reductions in 
coolant mass flow rates leading to localized thermal instabilities 

* Reductions in coolant mass flow rates result in further increases in local temperatures exacerbating changes to 
coolant fluid properties leading to localized thermal runaway 

* Unchecked localized thermal runaway leads to localized fuel melting 

Reactor designs with randomized flow paths are vulnerable to localized thermal instabilities, localized thermal runaway, 
and localized fuel melting 

Reactor designs and concepts can be assessed for thermal stability using the terms 

(=S) + (+S) = S T (=S)/V (+S)/V (S T )/V 

where 

=S Surface Area that is Parallel to the AP 

+S Surface Area that is Perpendicular to the AP 

S T Total Surface Area 

V Volume 

Some thermal instabilities can be eliminated with the use of cross flow pumps to establish a small perpendicular AP 
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